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Abstract 

We report exact analytical expressions locating the — > 1, 1 — > 2 and 2^4 bifur- 
cation curves for a prototypical system of two linearly coupled quadratic maps. Of 
interest is the precise location of the parameter sets where Naimark-Sacker bifur- 
cations occur, starting from a non-diagonal period-2 orbit. This result is the key to 
understand the onset of synchronization in networks of quadratic maps. 
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An essential clue to understand how dynamical systems synchronize and the 
origin of their riddled basins of attraction is provided by the mechanisms 
underlying the transverse desestabilization of low periodic orbits embedded 
in chaotic attractors[l,2,3,4]. So far, these mechanisms have been investigated 
only via direct numerical simulation of the dynamics. The purpose of this 
paper is to report exact analytical results for a prototypic model, two linearly 
coupled quadratic maps, defined by the following pair of equations 

xt+i = a - + b {xt - yt), yt+i = a - yj + h {yt - Xt). (1) 

Here Xt, yt represent dynamical variables, a is the local nonlinearity, h the 
coupling, and t = 0,l,2,... is the discrete time. Systems of coupled nonlinear 
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oscillators are nowadays used in practical applications, e.g., in the investiga- 
tion of coupled p-n junctions [5] and of Josephson-junction arrays [6]. They 
also serve as useful approximations of nonlinear partial differential equations 
describing, for example, the onset of fluid turbulence [7] and reaction diffusion 
systems [8]. Previous analytical results were obtained using linear stability 
analysis[9,10,ll], and renormalization thcory[12]. Here we go well beyond this 
presenting (i) the algebraic varieties delimiting stability domains for motions 
of periods 1 and 2, and, (ii) the Naimark-Sacker line (also called Hopf bi- 
furcation line) showing where quasiperiodicity stabilizes following a period-2 
oscillation. The motivation for considering this problem anew is our interest 
in using more realistic and flexible maps as local oscillators in applications 
dealing with diffusive-advective effects like, e.g. in ocean convection[13,14,15]. 
For additional discussion of problems involving coupled systems we refer the 
reader to the excellent tutorials in Refs. [1,2,3,4]. 

As is well-known [16, 17], there are three ways in which a flxed point p of a 
discrete map / may fail to be hyperbolic: Df{p) may have an eigenvalue -|-1, 
an eigenvalue —1, or a pair of complex eigenvalues, say A and A, with |A| = 1. 
The flrst two situations are familiar from period-doubling cascades and appear 
much more frequently than the last one. But it is precisely this last situation, 
characteristic of a Naimark-Sacker bifurcation, that is of great interest for 
coupled systems and that we consider here. 

Equation (1) has four flxed points, two diagonal and two non- diagonal: 

Di = {xi,Xi), D2 = {X2,X2), Ni = {x3,Xi), N2 = {x4,Xs), (2) 

where xi = (-1+ 0/2, X2 = (-1-0/2, = {2b-l+r])/2, x^ = {2b-l-r])/2, 
and = Vl + 4a and r] = Vl + 4a — 46^. The stability of these points is ruled 
by the equation | J — A/| = 0, where A is the eigenvalue, I is the identity matrix 
and J is the Jacobian matrix of the mapping. The eigenvalues of these four 
flxed points are always real, with D2, Ni and N2 being always unstable since 
at least one of their eigenvalues is larger than 1. The remaining flxed point, 
Di, is stable in the region delimited by —1/4 < a < 3/4 and \/l + Aa/2 — 1 < 
h < \/l + 4a/2, bounded by the curves 

1 3 

a = --, a=-, ci = 4a -46^ + 1 = 0, C2 = 4a - 46^ - 3 - 86 = 0, 
4 4 

shown in Fig. 1. 

The stable flxed point bifurcates into two quite different period-2 orbits when 
one increases a, depending on the value of b. When crossing the line a = 3/4 

one flnds diagonal period-2 orbits with signature {u, u) — ^ (f , v) {u, u) 
. . . , where n = (1 + V-3 + 4a) /2 and v = {I - V-3 + 4a)/2. We denote 
such orbits by P2D. Beyond the curve C2 = one flnds non diagonal or- 
bits with a signature (s, t) (t, s) — > (s, t) — > ... where s — (26-1-1-1- 
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Fig. 1. Stability domains for period-1 and period-2 orbits. P1D and P2 refer to do- 
mains of diagonal and non-diagonal orbits, respcctiA'oly. Nunil)(^rs indicate periods. 

V4a-462 -86-3)/2 and i = (26 + 1 - V4a - 46^ - 86 - 3)/2. These we 
represent by PI. 

The P2D orbit has real eigenvalues and is stable in the domain delimited by 
the curves 

3 5 

a = -, a = ^, C3 = 4a -46^ - 3 + 46 = 0, C5 = 4a - 46^ - 5 + 46 = 0. 

It loses stabihty along C5 = or a = 5/4, after a 2 — > 4 bifurcation. There is 
a P4L' diagonal orbit when crossing the line a = 5/4, and a P4 non-diagonal 
orbit when crossing C5 = 0. 

The eigenvalues of the P2 orbit are 



and can be complex. From these expressions we obtain[16] the following for- 
mula for the Naimark-Sacker bifurcation line: 



66^ + 106 - 4a + 4 ± 26^562 + 106 - 4a + 4 



(3) 




(4) 



It is important to observe that the linear coupling of the oscillators was fun- 
damental in allowing the derivation of the above analytical result. Already a 




Fig. 2. Stability domains obtained by discretizing the parameter interval in a mesh 
of 1200 X 600 points. Different gray shadings represent stability domains of attractors 
with different periods. 

quadratic coupling generates formulas which are intractable with present day 
hardware and commercially available software. 

From the eigenvalues along C2 = as well as the condition Ai = A2 = 1 we 
find that: (i) for b > 0, an unstable P2 orbit is born along C2 = stabilizing 
however along C4 = 4a — 46^ — 3 — 126 = 0, and (ii) for —1 < 6 < 0, a stable P2 
orbit is born, also along C2 = 0. Figure 1 summarizes all period-2 bifurcation 
loci. In the small diamond-shaped region in Fig. 1, delimited by C3 = 0, C4 = 0, 
C5 = and the Naimark-S acker line (indicated by NS), there is coexistence of 
P2 and P2D motions. On the NS line the P2 orbit disappears, to give rise to 
two closed curves (limit cycles). 

The stability domains for the coupled quadratic maps can be generated nu- 
merically, using the method described in Ref. [18], and they are shown in 
Fig. 2, displaying the richness of the several coexisting motions in parameter 
space. In this figure the basin of unbounded (diverging) attractors is indicated 
by 00. The lower border line seen between the P2 domain and black domain 
is defined by the Naimark-Sacker line. This means that inside the black do- 
main, immediately below the P2 region, we have quasiperiodic motion with 
two limit cycles originated from P2 trajectories. As we penetrate more and 
more into this black region, we reach a region of chaos. The same happens 
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Fig. 3. Quasiperiodicity and chaos seen in the phase-space of two hnearly coupled 
quadratic maps. For a = 0.5 and {xo,yo) = (0.1234,0.005), 1000 points have been 
plotted, after a transient of 400 iterations. Shown are three quasiperiodic solutions, 
for b = —0.35, —0.38, —0.40. Here, smaller b corresponds to larger radius. A fre- 
quency-locked solution with pcriod-6 appears at 6 = —0.41 (not shown in the figure). 
The complicated dotted structure is a chaotic orbit for b = —0.46. 

for P4 oscillations, the difference in this latter case being that the region of 
quasiperiodic motions has four limit cycles originated from P4 orbits. 

Figure 3 shows the birth and evolution of the Naimark-Sacker bifurcation in 
phase space when we walk along the vertical line a — 0.5 in parameter space. 
Here we plotted the amplitudes of the maps for four representatives values 
of the parameter b. For b close to the Naimark-Sacker line the motion covers 
a closed curve and is quasiperiodic. For smaller values of b, a quasiperiodic 
transition to chaos takes place. Periodic windows, not visible in Fig. 3, appear 
in this route to chaos via quasiperiodicity. 

We demonstrated that two linearly coupled quadratic maps can show quasiperi- 
odic motion arising from a Naimark-Sacker bifurcation of a P2 orbit. This 
was shown analytically and corroborated through a numerical simulation, as 
may be seen from the similarity between Figs. 1 and 2. The first figure was 
obtained analytically while the second numerically. The precise analytical de- 
termination of the Naimark-Sacker road to stable quasiperiodic motions opens 
the possibility of several numerical experiments, in particular, to applications 
connected with synchronization in networks of maps[13,14,15]. Another inter- 
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esting application is to investigate how a mismatch in the parameters would 
affect the bifurcation structure. These applications will be reported elsewhere. 
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